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Abstract 

Disturbing of a spacetime geometry may result in the appearance of an 
oscillating and damped radiation - the so-called quasinormal modes. Their 
periods of oscillations and damping coefficients carry unique information 
about the mass and the angular momentum, that would allow one to 
identify the source of the gravitational field. In this talk we present recent 
bounds on the diffused energy, applicable to the Schwarzschild spacetime, 
that give also rough estimates of the energy of excited quasinormal modes. 
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1 Introduction 



Quasinormal modes are scattering-type solutions of the Schrodinger equation 
which satisfy a quite peculiar boundary condition - that at both "ends" the 
waves are purely outgoing. They are being studied in the context of general 
relativity (the Schrodinger equation emerges there through the standard sepa- 
ration of the time dependence in a wave equation) in the course of the last thirty 
years (f|]-||). Much is known about their eigenvalues and their temporal evolu- 
tion in the case of Schwarzschild black holes and neutron stars. An exhaustive 
review on that topic is H ; see also Q . An observer located at a fixed space po- 
sition would find that quasinormal modes oscillate in time and their amplitude 
exponentially decays. The period of oscillations and damping coefficients carry 
unique information about the mass and the angular momentum. 

Perturbation of a spherically symmetric spacetime geometry Q can be de- 
scribed - as far as the backreaction can be neglected - by a linear wave equation 
(|| - H ). The interesting fact is that a (compact support) perturbation of 
a spacetime geometry may result in the appearance of an outgoing radiation 
that coincides in a bounded region of spacetime with a linear combination of 
quasinormal modes. At some intermediate time - later on the so-called tail term 
dominates - the perturbation is dominated by the fundamental mode, since the 
latter is damped weaker than the other quasinormal modes. The oscillation 
periods and damping coefficients do not depend on perturbations; this feature 
can be used in order to identify the source of the gravitational field [Q . 

It is of interest to know how much energy can be carried by quasinormal 
modes. Their energy can be estimated by the so-called diffused energy - the 
energy loss that is due to the backscattering. In this paper we present recent 
results in this direction. The order of the rest of this paper is as follows: 

i) Space-time curvature and two patterns of propagation of massless fields; 

ii) Vibrations of a spacetime - an example; 

iii) Recent results on the energy diffusion in the Schwarzschild spacetime; 

iv) Dependence of backscatter on the frequency of waves; 

v) Discussion; 

vi) Lessons from numerics. 



2 What is backscattering? 

Let an outgoing null cone T a originate from (a,0). Assume that a flash of 
radiation is initially purely outgoing JlO| and that its support is contained in an 
annular region (a, 6), b < oo. Then, depending on whether or not the spacetime 
is curved, the following can be observed. 

i) In the flat Minkowski spacetime a wave remains purely outgoing. No radia- 
tion can be found in the interior of the cone (Fig. 1). In this case no backscatter 
(@i @) occurs (in the Hadamard's terminology: the type B Huyghens principle 
holds true). 
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Figure 1: Solid line represents an outgoing null cone. Arrows show the direction 
of the radiation. 

ii) A wave backscatters in a curved, e.g. Schwarzschild, spacetime. Some 
energy, denoted later as SE a , diffuses inward through T a and is lost from the 
main stream (Fig 2). 
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Figure 2: Solid curve represents an outgoing null cone. Arrows show the direc- 
tion of the radiation - there appears an ingoing component. 

One can shortly state that waves propagate in a curved spacetime like elec- 
tromagnetic waves in a medium with a varying refraction index. A fraction of 
the radiation scatters off the curvature of the geometry and a part of the initial 
energy never reaches infinity. 

3 Damped oscillations 

Figure 3 shows a picture that is typical in the case of perturbations that have 
compact support. An initial perturbation is depicted in Fig. 3(a). Some time 
after the initial pulse runs through the observer, he (or she) can observe that an 
oscillatory (single frequency) radiation dominates. As pointed out before, the 
period and damping coefficients are independent on the perturbation. 
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There are two interconnected problems that can made difficult the identi- 
fication of the dominant mode. First, in the asymptotic zone, t >> 2m, the 
tail terms dominate since they decay as some power of 1/t. Second, the 
amplitude of the oscillations quickly decreases, as exemplified by Fig. 3(b). 



(a) y/(t) 




Figure 3: Time evolution of the Regge- Wheeler function for L = 2. The 'ob- 
server' is located at r* = 800 * m 

4 Quantitative estimates of the backscattered 
energy. 

We specialize, beginning from this section, to spherically symmetric spacetimes. 
The backscattered (diffused) energy will be found from the energy conservation 
and some potential estimates. Let us point out that our intention is to derive a 
bound from above on the fraction of the backscattered energy; this bound should 
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be independent of the details of the initial pulse of the outgoing radiation. The 
only information that is required, is the mass m of the central object and the 
initial position a of the inner boundary of the outgoing radiation. (In the case of 
nonspherical spacetimes one presumably would need also the information about 
the total angular momentum.) 

The idea that it is possible to bound the backscattered energy entirely in 
terms of initial energy and initial position was first put on trial in the example 
of the massless scalar field (Q - [|lj|; a recent work is Q). The Theorem 
formulated below summarizes results that have been obtained while implement- 
ing a programme that was formulated in the case of electromagnetic radiation 
in flr| ) and then modified in [jlTj and [jl8|. While the analytic part con- 

cerning the electromagnetic radiation is fairly satisfactory (any further progress 
would most likely involve advanced numerical methods), the work on the grav- 
itational radiation is not finished. Preliminary results on the polar modes are 
given in |l7j while the investigation of axial modes is in progress. 

Theorem. 

Assume a flash of dipole electromagnetic radiation (quadrupole gravitational 
radiation - axial or polar) that is initially purely outgoing and that has support 
in an annular region (a, 6), b < oo. Let the initial energy be E a (0) and m = m/a. 

Then the fraction of the diffused energy SE a / E a (0) satisfies the inequality 



< C(m) - 



E a (0) ~ \ a J 



Here 

i) (the case of electromagnetic radiation) for a > 10y , 0.3 = C(0) < C(m) < 
(7(10/3) w 1.7 

Cirri) = 7— — r^— r- 1 -2760m 5 + 828m 4 + 

V ; 10080((-1 + 2m) 4 m 4 ) ^ 

44m 3 + 2352m 6 + 2016m 4 ln(l - 2m) + 
2688m 6 ln(l - 2m) - 4032m 5 ln(l - 2m) - 
360m 3 ln(l - 2m) + 6m + 36m 2 ln(l - 2m) - 

30m 2 + 3 ln(l - 2m) - 18m ln(l - 2m) J ; 

(1) 

ii) (axial gravitational waves ) JT9| 1.05 = C(0) < C(m) =? ; 

iii) (polar gravitational waves) Q 55/2 = C(0) < C{rh) =? (conjecture: 
C(0) < 4). 



Sketch of the proof 
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4.1 Preliminaries. 



A Schwarzschildean line element reads (we neglect the backreaction effect): 

2m 1 

ds 2 = -(1 - — )dt 2 + 2^dR 2 + R 2 dn 2 . (2) 

R 1 -~ 

The radial terms of the multipole expansion satisfy the following (reduced) wave 
equations (§,[§, §J), 

(-d 2 +d 2 ^ t =V(R)*t. (3) 
In the the case of lowest multipoles we have 

2m \ 2 
R >W 



i) I = 1 and V(R) = (1 — (dipole term, electromagnetism); 



ii) I = 2 and = (1 - ^ff ) ^1 - if J (quadrupole term, axial gravita- 

tional waves); 

hi) I = 2 and = (1 - ^1 - ^ + ^^^l^^ (quadrupole 

term, polar gravitational waves). 

In what follows we shall deal only with the electromagnetic case; the other 
two cases can be treated similarly. Although the calculations become more 
complex, basic scheme is the same. 

The most general solution of the dipole electromagnetic radiation in the 
Minkowski spacetime is given by 



b ( J(fi - t) - 9(fi+t) ) + zgza±£g±j> , (4) 

where the /-related part describes the outgoing radiation while g is responsible 
for the ingoing wave. For any initial data, one can uniquely determine / and g 
and in this way specify the outgoing and ingoing initial pulses. 

We invoke to the above decomposition in order to specify what is meant by 
in- or out- directed waves in the Schwarzschild spacetime. Define the Regge- 
Wheeler variable r*(R) = R + 2mln(^ — 1). Then having initial data, one can 
construct functions / and g in a similar way as in the Minkowski spacetime. If 
from the |l9| construction follows that 3 = 0, then we will say that the radiation 
is purely outgoing. In such a case it is useful to define 

*(R,t) = -a r .f(r*-t)+ f ^-*\ (5) 

and to seek a solution of the dipole wave equation, ^f(r*,t), having the following 
form 

* = * + 5. (6) 

We would like to point out that at t — the function 6 vanishes up to the first 
time derivative JIB], (5 = and OqS = 0. One finds that the evolution equation 



G 



reads 



(-5 o 2 + r 2 ,)* = (l-^) 



R 2 



6m/ 
R 4 



(7) 



4.2 Main steps of the proof. 



The energy of the electromagnetic field that is contained in the annulus (R(t) , oo) 

, ,., 2m, /a , T . 2 , 2(*) 2 \ 



4.2.1 Bound on / in terms of the initial energy. 

Lemma 1. Define nji = 1 — m = m/a and y = R/a. Let f(a) = 0, and 
a > 10m/3. Then 

J 2 (R,0) E a (0) i(mw( - , ... 

< — arj {R)F(m,y), (8) 



where 



i? 2 1 " 8tt 



16m 4 16m 4 

F{m,y) =y-l + 



3(-y + 2m) 3 3(-l + 2m) 3 
16m 3 16m 4 24m 2 



(-y + 2m) 2 (-l + 2m) 2 (-y + 2m) 

24m 2 y — 2m 

1 8m In ■ 



(-l + 2m) l-2m 



(9) 



This is essentially a Sobolev-type estimate. For the proof see [ fl8[ . 
4.2.2 Estimate of an "energy" of S 

This "energy" denoted as H is not conserved - there is a volume-dependent term 
in the integral form of this "energy" evolution law. It appears to be, however, 
a useful quantity. H is defined by 

H(R,t)= [ X 'dr( { -^+ Vr (d r 6) 2 + %); (10) 
Jr. v Vr r 2 ) 

this integral is done on a fixed Cauchy hypersurface t — const. One can prove 
Lemma 2. Let the support of initial data be (a, b), a > 10m/3, b < oo. 
Then 
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H(a t ,t) < 36ra 2 



ds 



dr 



1/2 



(11) 



where the i-integration follows along T a ^ at t ^ while the r— integration is done on 
a fixed Cauchy slice. 

The crucial point in the proof of Lemma 2 is that 



H(a t ,t) = - 



SE a {t) 
2tt 



- 12 / ds drd S 



f 



(12) 



dropping out the nonpositive 8E a terms and using the Schwarz inequality yields 
the lemma Jl8[ . 

4.2.3 The energy conservation and bounding of the diffused energy 
(the energy "loss") 

Define r^ (^. t ) - a segment of an outgoing null geodesic that connects (Ro,t = 
0) with (R,t). Further, let us introduce 



h-(R,t) = 



i 2m 

1 R 



(d + d r *)6; 



(13) 



this corresponds to this component of the "reduced" strength field tensor that 
is directed inward. 



The rate of the energy change along r a is given by 
(d Q +d r *)E a = 

2m,/ do* 



„ 2m. 

- 2 " (1 - -r' 



"-f>(r 



2m 
R 



-2tt(1 



2m 
~R~' 



(1- 



2m, 



h_ - 



R 2 



2 



(14) 



but / = * = while their first derivatives vanish on T a - The integration of 
4j) along f gives the energy that diffused through the segment T a .a{t), 



SE a (t) = E a (Q) - E a (t) 

ca{t) 



2TT 



dr 



2m, 2S 2 
(1 - — ) h - + — 



^From the ^-conservation law and Lemma 2 one gets 



SEa 

2vr 



< 36™ 



ds\ 



dr 



/V 



l/2-i 2 



(15) 



(16) 



here SE a = lim t _ >00 6E a (t). 



The " electromagnetic" part of the main Theorem follows from the preceding 
bounds. 



5 Dependence of the backscatter on the frequency 
of waves 

The backscattering depends on the relative width of support. That is well 
known from the numerical analysis of Vishveshwara pl]j), but the first - up to 
our knowledge - proof appeared in |L6| . Below we sketch the main results. When 
support of the initial radiation is very narrow, k = (b — a)/a << 1, then 



6E a <c (2m\\ 

In the limit k — > ^{o) ~ * ^> * ne backreaction is negligible. We would like to 
point that this argument works for all a > 2m, in contrast to the main Theorem. 

The question arises: how do we interpret this? The similarity theorem p2| 
of Fourier transform theory states that compression of support of a function 
corresponds to expansion of the frequency scale. It means that 

,. E(w < fU 

lim v ~ r ; -> 0: 17 
E a (0) 1 V ' 

above E(u> < fl c ) is the electromagnetic energy of modes with frequencies 
smaller than a fixed frequency, lj < VL C . Thus we can conclude that the high 
frequency radiation is essentially unhindered by the backscattering while long 
waves can be backscattcrcd. 



6 Discussion 

The diffused energy 5E a bounds from above the sum SE qmt +SEf where 5E qmt 
is the energy carried by the quasinormal modes and the tail, and SEf is the 
energy of that radiation that falls to a black hole or hits the surface of a star. 
Below we present several data that estimate from above the total fraction of 
the backscattered energy. A more detailed calculation would show that the last 
contribution, 5Ef, dominates, so that we expect that SE qm t/E a (0) can be a 
small fraction of the number given in forthcoming examples. 

In the first group of examples we consider the case, when the initial pulse is 
close to the Schwarzschild radius R = 2m; 
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i) a = 10to/3 (e. g., a surface of a supercompact neutron star ): e (o) ^ ^-5; 

ii) a — Am (e. g., a supercompact neutron star ): e E (o) ^ 0-3; 

iii) a = 5m (a standard neutron star): -J^fjy < 0.13. 

It is instructive to notice that this effect is very weak in the case of other 
astronomical objects. For the Sun for instance : gPfo) ~ 10 -13 , while for white 

dwarves: e B (o) < 10~ 7 - Our estimates allow one also to answer the question 
that was raised in the literature (see, e.g., Q), how strongly gravitational field 
of a galactic interior can damp the outgoing radiation? One can infer from the 
preceding examples that the effect may matter only, if there is a black hole in 
the interior of a galaxy. 



7 What can be learned from numerics? 

We summarize shortly the main conclusions that can be obtained through a 
numerical analysis. 

i) If a radiating source is close to a horizon, then the damping can be quite 
strong. We found, in accordance with expectations, that the backscatter is 
strong deep inside the photon sphere. There are known examples in which from 
20 percent (electromagnetism) to 49 percent (gravitational polar modes) of the 
initial energy gets diffused |jl8); that should go to 100 % if a very compact source 
is infinitely close to the event horizon. 

ii) The efect depends on the width of a"renormalizcd" ((b — a)/r] a ) support 
of initial energy. The backscatter is strongest, when the width is of the order of 
a - comparable to the areal distance from the center. 

iii) The present analytic estimates are not exact. They are expected to yield, 
in combination with appropriate numerical methods, more satisfactory results. 
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